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VARIATIONS OF THE LUNAR ORBITAL PARAMETERS 
OF THE APOLLO CSM-MODULE 

by 

T. P. Timer 

Goddard Space Flight Center 


SUMMARY 

This report investigates the variation of the orbital parameters for the 
forty hour period after insertion into two reference Apollo Lunar Parking Or- 
bits* . It is found that while the variation of most of the parameters under the 
influence of a triaxial Moon and other perturbing forces is small, the change in 
pericynthion height is significant. In Table 1 the amount of change in the various 
elements during this time is shown. 

The most pertinent piece of information derived from this investigation is 
with regard to the choice of hypothetical lunar landing sites on the surface of the 
Moon (see Figure 1). Choosing a point on the earthward side of the trace which 
satisfies the constraint that it has a longitude of between 315° and 45° and a lati- 
tude of +5° and -5° with respect to a selenographic coordinate system one has a 
hypothetical landing site. It is found that any subsequent trace of the sub- module 
point at the same longitude will always lie within 0°5 in latitude of the initial 
landing site chosen during the forty hours under consideration. Therefore, 
landing sites chosen in the above manner meet the two requirements as de- 
scribed in Apollo Navigation Technical Report 2 which states that the two con- 
straints on the lunar parking orbit mission are as follows: 

(1) the plane of the trajectory must contain a landing site vector on the 
earth side of the Moon which has a longitude of between 315° and 45° 
and a latitude of between +5° and -5° in selenographic coordinates. 

(2) the plane of the orbit must be oriented so that the lunar landing site 
doesn't move out of the orbital plane more than 0°5 during the period 
of 3 to 39 hours after lunar insertion. 



*pp. 3*10 Apollo Navigation Group Technical Report 65*AN-1.0, February 5, 1965. 
2 pp. 3-6, 3-7. 


ill 


VARIATIONS OF THE LUNAR ORBITAL PARAMETERS 
OF THE APOLLO CSM-MODULE 


by 

T. P. Timer 

Goddard Space Flight Center 


INTRODUCTION 


The basic concept behind this report was to see how the orbital parameters 
for the two reference lunar parking orbits behaved during the forty hour period 
after lunar insertion. To this end two 7090/7094 computer programs were used 
to generate this information which is presented as two series of graphs. Series 
A graphs are associated with the 1968 reference trajectory. Series B graphs 
are associated with the 1969 reference trajectory. It should be noted at the out- 
set that several coordinate systems are in use throughout this report. The 
three coordinate systems in use are the following: 

center: selenocentric 

fundamental plane: parallel to the earth’s equatorial plane 
x-axis: mean equinox of date 

center: selenocentric 

fundamental plane: true equator (earth) of date 
x-axis: true equinox of date 

center: selenocentric 

fundamental plane: true equator of the moon 

x-axis: intersection of the prime meridian and the 
equator of the moon. 

The two computer programs used as sources of data are the Jet Propulsion 
Laboratory's Space Trajectory program and the Goddard Interplanetary Tra- 
jectory Encke Method (Version 12-C) program. Both programs numerically 
integrate the orbit by the Encke Method hence one program serves as an inde- 
pendent check of the validity of the data generated by the other program. 
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I. General Description of the Reference Trajectories 


The lunar parking orbit is constrained to be near- circular with a semi- 
major axis of 1018.5 nautical miles (1886 km) with an altitude above the lunar 
surface of 80 nautical miles ±5 nautical miles (148 ± 9 km). The duration of 
this trajectory is constrained to be between two to sixty hours 3. In this study 
both cases meet these requirements. The insertion conditions for the two 
reference 3 4 trajectories under consideration (study) are: 


Date of Lunar Insertion Date of Lunar Insertion 

Reference Trajectory #1 Reference Trajectory #2 


May 9, 1968 

15 hrs/54.0 mins/56.610 secs. 


September 20, 1969 
5 hrs/10 mins/12.176 secs. 


X 

= 

-1012.437 

0 

n.mi. 

X 

=5 

+ 165.639 

36 n.mi. 

Y 

= 

- 51.263 

872 

n.mi. 

Y 

= 

- 919.376 

95 n.mi. 

Z 


- 20.120 

039 

n.mi. 

Z 

= 

- 415.861 

51 n.mi. 

X 


- 287.896 

06 

ft/sec 

X 

— 

-5208.366 

5 ft/ sec 

* 

= 

+4914.267 

5 

ft/ sec 

Y 

= 

- 817.043 

16 ft/sec 

z 

= 

+1967.337 

7 

ft/sec 

Z 

“ 

- 268.213 

99 ft/sec 

X 


-1875.033 

324 

km 

X 

— 

+ 306.764 

095 km 

Y 

= 

- 94.940 

691 

km 

Y 

= 

-1702.686 

111 km 

z 

= 

- 37.262 

312 

km 

Z 

= 

- 770.175 

517 km 

X 

— 

- 0.087 

750 

895 km/sec 

X 

- 

1.587 

513 284 km/ sec 

Y 


+ 1.497 

871 

730 km/sec 

Y 

= 

0.249 

035 253 km/ sec 

Z 

= 

+ 0.599 

645 

730 km/sec 

Z 

= 

0.081 

751 788 km/sec 


conversion factors: p 4-8 


COORDINATE SYSTEM 


center: 
fundamental plane: 
x-axis: 


selenocentric 

plane parallel to earth's mean equator of date 
mean equinox of date at insertion. 


3 P 3-7. 

4 p 3-11. 
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In terms of orbital elements, the insertion conditions, using the same 
coordinate system above, are: 

Reference Trajectory #A Reference Trajectory #B 

1968 1969 



a = 

1878.569 

2 

km 

a = 

1894.578 

3 

km 


e = 

0.000 

406 

737 41 

e = 

0.000 

419 

706 10 


i = 

158°183 

342 


i = 

155°804 

726 


« 

n = 

180!056 

736 


a = 

182?414 

087 



CO = 

357J743 

31 


CO = 

262°877 

900 



m = 

359H.96 

573 


m = 

0?000 

824 

8 


EL. Discussion of the Graphical Results 

The series of graphs A-l through A- 6, and B-l through B-6 are related to 
the reference trajectories A and B respectively and show similar trends. Hence 
only the first set of graphs will be discussed in detail with much of the dis- 
cussion being applicable to both sets. 

By inspecting graphs A-l to A- 4 the graphs of the semi-major axis, eccen- 
tricity, inclination, and longitude of the ascending node it can be seen that the 
orbit chosen is relatively a very stable orbit. The semi-major axis varies at 
most about 0.5 kilometers from its original length during the 40 hours of in- 
vestigation. The inclination changes less than 0.05 degrees. The eccentricity 
is plotted to a scale with a factor 10“ 3 . Its fluctuations during the forty hour 
period are of the order of magnitude of 0.3 x 10~ 3 as a maximum. The longitude 
of the ascending node for all practical purposes changes linearly with respect 
to time, varying less than 0.05 degrees during the time interval considered. The 
purpose of the preceding discussion was to show quantitatively how these various 
elements varied with time. Figure A- 5 and Figure A- 6, however, provide the 
most useful information directly affecting a lunar parking orbit. In Figure A-5, 
height (radius vector minus the mean lunar radius) versus hours from epoch, 
shows how the height varies periodically with the period being the time for one 
circumlunar orbit. This variation is less than 0.15 kilometers. Figure A-6 
shows the limiting cases, the first and the last orbit data, of the trace of the 
command module on the lunar surface or the lunar sub-module point. It can be 
seen that for the worst case by choosing a particular point on the trace of the 
first trajectory and then choosing a point on the trace of the last trajectory at the 
same longitude, the variation in latitude would be at worst 0°2. A considerable 
portion of the time, this variation is less than 0!2. In graph B-5 a small secular 
difference is evident in the value of the ascending node as computed by the JPL 
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and Goddard Item programs. The reason for this slight difference has not been 
ascertained except that it is not due to numerical differences in the values of the 
physical parameters used by the program. (See Appendix A) 

CONCLUSIONS 

(1) Based on the graphical data presented in this report, either the 1968 or 
1969 reference trajectory could provide an orbit stable enough to meet the mis- 
sion requirements of having suitable landing sites on the lunar surface which do 
not vary more than 0?5 in latitude from the original landing site. This informa- 
tion is presented in graphs A- 6 and B-6 respectively where arbitrary hypothetical 
landing sites have been chosen. 

(2) The use of two computer programs, one to serve as an independent 
check on the data generated by the other, is justified as the inter comparison 
of the data on the various graphs indicates. 

(3) During this lunar parking orbit phase it would seem that the orbital 
path may be considered a constant with regard to three station USBS ground 
tracking. This statement is based on the knowledge that 3o- errors in position 
including both measurement noise and station bias errors are 8,500 meters. 5 
The change in the semi-major axis as found by this report is at most 500 meters 
which would be undiscemible since it is less than the accuracy of the position 
measurements made. 

In Table 1 the maximum changes in the parameters for the two cases are 
shown. Essentially, the results are identical. 


5Apollo Navigation Working Group Technical Report No. 65-AN-2.0, pp. 6-5, September 1 , 1965. 


Table 1— Maximum Changes in Orbital Parameters 



Aa 

Ae 

Ai 

Case #1 




1968 Data 

0.5 km 

0.3 x 10" 3 

0:05 

Case #2 




1969 Data 

0.45 km 
1 

0.3 x 10" 3 

0:05 



An 

AHp* 

A Latitude 

Case #1 
1968 Data 

0°05 

0.15 km 

0^2 

Case #2 
1969 Data 

0!20 

0.20 km 

O 

• 0 

cn 


*Hp “ pericynthion height. 
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Figure 1 
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Fiirure A- 1 Semi- major Axis (kilometers) Vs Hours from Epoch 
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APPENDIX A 


CONSTANTS 






ANWG 

ANWG 

JPL 

NAME 

QUANTITY 

JPL VALUE 

GSFC ITEM VALUE 

APOLLO 

NAME 

APOLLO VALUE 
& QUANTITY 

GME 

gravitational coefficient for 
earth 

3.986 0320 x 10 s km 3 /sec 2 

19.909 4165 er 3 /hr 3 

n EARTH 

gravitational parameter for 
earth 

3.986 032 x 10» 4 mVsec 2 
1.990 941 65 x 10 eH/hr 2 

GMM 

gravitational coefficient for 
moon 

4.900 7589 x 10* kra*/sec 2 

0.244 782 8913 erVhr 2 

n MOON 

gravitational parameter for 
moon 

4.902 778 x 10 12 mVsec 2 
0.244 883 757 erVhr 2 

GMS 

gravitational coefficient for 
sun 

1.327 154 4 x lo 11 krn^sec 2 

0.513 736 472 x 1 Q-* auVhr 2 

u SUN 

gravitational parameters for 
S' 111 

1.327 154 45 x 1 O 20 m 3 /sec 2 
6.628 865 68 x erVhr 2 

GMA 

GMV 

gravitational coefficient for 
mars 

gravitational coefficient for 
venus 

4.297 779 9 x lO^kmVsec 2 
3.247 695 0 x 10 5 kmVsec 2 

0.166 365 411 7 x lO -12 au 3 /hr 2 
0.125 717 034 3 x 10' u au 3 /hr 2 

n.a. 

n.a. 

n.a. 

n.a. 

GMB 

gravitational coefficient for 
bar y center 

4.035 039 5 x 10 5 km 3 /sec 2 

not in GSFC Item 

n.a. 

n.a. 

GMJ 

gravitational coefficient for 
jupiter 

1.267 106 0 x 10 8 kmysec 2 

0.490 492 590 7 x l (T 9 au*/hr 2 

n.a. 

n.a. 

J 

coefficient of the second 
harmonic in earth’s oblat- 
ness 

1.623 4500 x 1(T 3 
non-dimensional 

1082.3 x 10-6 
non- dimensional 

J 

C 2.0 

constants associated with the 
gravitational potential of the 
earth 

1.623 45 x 10' 3 
-1082.30 x 10" 6 

H 

coefficient of the third 
harmonic in earth's oblat- 
ness 

-0.574 999 99 x 10* 5 
non-dimensional 

-2.3 X 10 r* 
non- dim ensional 

H 

t-3 . o 

constants associated with the 
gravitational potential for the 
earth 

-0.575 x 10‘ 5 
2.3 x 10 -6 

D 

coefficient of the fourth 
harmonic in earth’s obial- 
ness 

0.787 499 99 x 10~ 5 
non-dimensional 

-1.8 x IQ - * 
non-dimensional 

D 

C 4.0 

constants associated with the 
gravitational potential of the 
earth 

0.787 5 x 10* 5 
1.8 x 10" 6 

RE 

earth’s radius to be used 
in earth's oblate potential 

6378.165 km 

6378.165 km 

Re 

earth radius 
6378.165 km 

G 

universal gravitational constant 
for lunar oblatness 

4.902 643 x 10 12 m 3 /sec 2 

0.244 782 9 er 3 /hr 2 

,* MOON 
mMOON 

gravitational parameters for 
the moon 

2.448 837 57 x 10* 1 er 3 /hr 2 
4.902 778 x 10 12 m 3 /sec 2 

A 

moment of inertia for the moon 
to be used in the lunar oblate 
potential 

0.887 459 98 x 10 29 kg-km 2 

0.887 46 x 10 35 kg-m 2 

A 

moment of inertia about 
principal rotational axis 
0.887 817 983 4 x 10 35 kg-m 2 

B 

moment of inertia for the moon 
to be used in the lunar oblate 
potential 

0.887 63 9 98 x 10 29 kg-km 2 

0.887 64 x 10* s kg-m 2 

B 

moment of inertia about 
principal axis 

0.888 001 9542 x 10 35 kg-m 2 

C 

moment of inertia for the moon 
to be used in the lunar oblate 
potential 

0.888 0099 8 x 10 29 kg-km 2 

0.880 1 x 10 35 kg-m 2 

C 

moment of inertia about 
principal axis 

0.888 369 78 1 7 x lO^kg-m 2 

OME 

rotational rate of the earth 

0.004 178 074 1 deg/sec 

0.262 516 170 81 rads/hr 

w 

angular velocity of the earth’s 
rotation 

4.178 074 16 x 10' 3 deg/sec 
0.262 516 142 rad/hr 

au 

astronomical unit to convert 
planetary ephemerides to km 

1.495 990 0 x 10® km/au 

1.495 99 x 10 8 km/au 

au 

astronomical unit ,, 
1.49 5 990 00 X 10" m 

REM 

earth radius to convert lunar 
ephemerides to kilometers 

6378.165 km 

6378.165 km 

n.a. 

n.a. 


b.i. = ikx applicable. 




